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Abstract 

An analysis is given of the decay fj, — > e + 7 in an MSSM extension with a vectorlikc generation. Here 
mixing with the mirrors allows the possibility of this decay. The analysis is done at the one loop level with 
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1 Introduction 


Lepton flavor violation provides a new window for physics beyond the standard model. Since there is no 
CKM type matrix in the charged leptonic sector, flavor violations involving charged leptons arise via loop 
corrections which in particular can produce charged lepton flavor violating processes such as —» £j/- 
Recently the MEG experiment [T] has put the most stringent bound thus far on the lepton flavor violating 
decay [i —> ey so that 

B(n -» e + 7) < 2.4 x 1CT 12 at 90% CL (MEG). (1) 

In this work we explore the implications of a new leptonic vector generation for the —> e + 7 decay. Specif¬ 
ically we consider an additional generation of leptons and their mirrors that mix with the three ordinary 
generations of leptons. Inclusion of a vectorlike generation brings in new sources of CP violation which enter 
in /r —> ey decay. These arise from diagrams where one has charginos and sneutrinos and neutralinos and 
charged sleptons in the loops. Additionally one has diagrams with W and neutrinos and Z and charged 
leptons in the loops. Such diagrams can produce observable effects and thus the experimental upper limit 
constrains the parameter space of models. Specifically we will show that the /i —> ey process can allow one 
to probe new physics arising from the MSSM extension. The reason for considering a vectorlike leptonic 
generation is the following: First vectorlike generations naturally appear in a variety of grand unified models, 
string models and D brane models and some of these can survive down to low scales |Sj- Second a vector like 
generation is anomaly free so the good properties of the model as a quantum field theory are protected. 
In previous works [3 si sum Mas ni muni we have considered the effects of an extra vectorlike 
generation on a number of processes and here we extend the analysis to discuss /r —► ey decay which is one 
of most stringently constrained lepton flavor violating process. We also investigate the effect of CP phases 
on the decay ji —>• ey. Vectorlike multiplets have also been used by other authors (see, e.g., [13 HSl USD- 
Further, /r —> ey decay has been analyzed in several previous works (see, e.g., mmmmmmmy 
However, none of the previous works explore the class of models discussed here. 

The outline of the rest of the paper is as follows: In section 2 we define the model with an extra vectorlike 
leptonic generation and specify the nature of mixings between the extra vectorlike generation and the three 
ordinary generations. In section 3 we give the interactions of the leptons and mirror leptons with the 
charginos and the neutralinos in the mass diagonal basis. In section 4 we give an analysis of the interactions 
of the leptons and their mirrors with the W and Z bosons. An analytic analysis of fi —> ey decay is given in 
section 5 which includes charginos and neutralinos in the loops as well as W and Z bosons in the loops. A 
numerical analysis of the /r —> ey branching ratio is given in section 6. Here it is shown that the vector like 
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generation gives a significantly large contribution which allows one to probe and constrain the extended 
model. It is known that CP phases can have a large effect on SUSY loop corrections (for a review see [24] ) 
and thus the effect of CP phases on the decay /x —► eq is also analyzed. Conclusions are given in section 7. 
Details of the scalar mass squared matrices are given in section 8. 

2 Extension of MSSM with a Vector Multiplet 

In this section we extend MSSM to include a vectorlike generation which consists of an ordinary fourth 
generation of leptons, quarks and their mirrors. As mentioned in section 1 vectorlike multiplets arise in 
a variety of unified models some of which could be low lying. In the analysis below we will assume an 
extended MSSM with just one vector multiplet. Before proceeding further we define the notation and give 
a very brief description of the extended model and a more detailed description can be found in the previous 
works mentioned above. Thus the extended MSSM contains a vectorlike multiplet. To fix notation the three 
generations of leptons are denoted by 


tpiL = 





(1,1,1); (1,1,0); * = 1,2,3 


( 2 ) 


where the properties under SU(3)c x SU( 2)l x U(1)y are also exhibited. The last entry in the braces such 
as (1,2,—1/2) is the value of the hypercharge Y defined so that Q = T 3 + Y. These leptons have V — A 
interactions. We can now add a vectorlike multiplet where we have a fourth family of leptons with V — A 
interactions whose transformations can be gotten from Eq.(??) by letting i run from 1 to 4. A vectorlike 
lepton multiplet also has mirrors and so we consider these mirror leptons which have V + A interactions. 
The quantum numbers of the mirrors are given by 


x c = 


( 1 , 2 ’ 2 ) 5 


£7 l ~(1,1,-1); N l ~( 1,1,0). 


( 3 ) 


Interesting new physics arises when we allow mixings of the vectorlike generation with the three ordinary 
generations. Here we focus on the mixing of the mirrors in the vectorlike generation with the three genera¬ 
tions. Thus the superpotential of the model allowing for the mixings among the three ordinary generations 
and the vectorlike generation is given by 
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W = + CijlhHWi Ti + AH^lKl + hH\x cj N L + mix ci E L 

+ /ii + KHW^U + ^H{^ eL e c L + h' 2 Hlr eL v c eL ] 

+ h^r^L + fcijX c %L + hr c L E L + Uv c tL N l + UIE l + f^ L N L 
+ f'^ioX^L + f”e c L E L + f'^ c eL N L , 


( 4 ) 


where ~ implies superfields, ip R stands for ip 3 L, tp^L stands for ip 2 L and ip eR stands for ip\L- In ?? we have 
suppressed terms such as e^E^ etc for simplicity. Their inclusion will not change our analysis substantially. 
The mass terms for the neutrinos, mirror neutrinos, leptons and mirror leptons arise from the term 


C = — 


1 

2 


d 2 W 

dAidAj 


ipiipj + H.c., 


( 5 ) 


where ip and A stand for generic two-component fermion and scalar fields. After spontaneous breaking of 
the electroweak symmetry, ((Hi) = V \/\/2 and ( H 2 ) = we have the following set of mass terms 

written in the 4-component spinor notation so that 


~ £ m =iR(M f )t L + riR(Me)ri L + H.c., (6) 

where the basis vectors in which the mass matrix is written is given by 

= (fo-fl N r v^ r v eR ) , 

= (y rL N l v eL ) , 

Vr = (tr Er /Ir e R ) , 

vl = (■ t l E l h l e L ) , (7) 

and the mass matrix Mf is given by 


M f = 


(f[v 2 /V2 

~h 

0 

V o 


h 

f 2 v i/\/2 

/s 

n 


0 

Kvi/y/2 

0 


0 \ 

-fs 
0 

h' 2 v 2 /V2) 


We define the matrix element (22) of the mass matrix as TOat so that 


m N = f 2 v i/a/2 - 


( 8 ) 


( 9 ) 


The mass matrix is not hermitian and thus one needs bi-unitary transformations to diagonalize it. We define 
the bi-unitary transformation so that 
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D R (Mf)D v L = diag (rn^, 1 ,m^ 2 ,m^ 3 ,m^ i ). 
Under the bi-unitary transformations the basis vectors transform so that 


( 10 ) 


( VTR \ 


/Vl B \ 

Nr 

U VR 

itti 

Q 

II 

^3 R 

Ve R J 


\V>4 R J 


^ v t l \ 

i 


N l 

"hr 

II 

V>3 x. 

V'aJ 




( 11 ) 


In ipi,ip 2 ,ip 3 iip 4 are the mass eigenstates for the neutrinos, where in the limit of no mixing we identify ipi 
as the light tau neutrino, ip 2 as the heavier mass eigenstate, ip z as the muon neutrino and ip 4 as the electron 
neutrino. A similar analysis goes to the lepton mass matrix Mp where 


M e 


(fiVi/y/2 

h 

0 

V 0 


In general / 3 , / 4 , f 5 , ft, ft, ft, ft, ft, ft can 


h 

0 

0 ^ 

ftv 2 /V 2 

ft 

ft' 

ft 

hiVi/y/2 

0 

ft 

0 

h 2 v 1 /V 2 ) 


be complex and we define their phases so that 


( 12 ) 


fk = |,/fc|e* Xfc , f k = \f k \e ix ' k , /"=|/"|e^' ;* = 3,4,5. (13) 

We introduce now the mass parameter niE defined by the (22) element of the mass matrix above so that 

m E = ftv 2 /V2. (14) 

Next we consider the mixing of the charged sleptons and the charged mirror sleptons. The mass squared 
matrix of the slepton - mirror slepton comes from three sources: the F term, the D term of the potential and 
the soft SUSY breaking terms. Using the superpotential of Eq. (4), the mass terms arising from it after the 
breaking of the electroweak symmetry are given by the Lagrangian 

C = Cp + Cp> + C so ft , (15) 

where Cf is deduced from Fp = dW/dAp and — Cf = Vp = F t F* i s given in the appendix, while the Cp is 
given by 

-C D = | cos 2 9 W cos 2^{v tL v* l - t l tI + + ^lKl - e l&l 

+ ErE* r - N r N* r } + * ro| sin 2 9 W cos 2P{v tL v* l + t l t* l + Vp L ft L + AlAI 
+ u eL v* eL + e L e* L - E r E* r - N R N* R + 2 E L E* L - 2f R f* R - 2 fi R ft R - 2e R e* R }. (16) 
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For £ so f t we assume the following form 


= M? L r T *Ji L + M\x ci *x ci + M? L ^; L ^ L + Mt L r e *Ji L + M?y T * L ^ TL + MlvZ L i>° L 

+ MIvTlKl + M?f c L *f c L + + M?eZel + M%E* L E L + M%N* L N L 

+ eMArHl^fl - f[A u Mrt L K L + hiA^Hl^l - KA^H^ l ^ l 
+ h, 2 A e Hlft eL e c L - h' 2 A v H$i L D c eL + f 2 A N H\x c °N L - f'A E H^x CJ E L + H.c.} . (17) 

Here Mg^, etc are the soft masses and A e ,A Ve etc are the trilinear couplings. The trilinear couplings 
are complex and we define their phases so that 

H e = \A e fe iaAe , A Ue = \A v<s \e iaA *e , ■ ■ ■ . (18) 

From these terms we construct the scalar mass square matrices. These are exhibited in section 8. 

3 Interactions with charginos and neutralinos 

In this section we discuss the interactions in the mass diagonal basis involving charged leptons, sneutrinos 
and charginos. Thus we have 

2 8 

~£r-p- x - = E E ^{C^Pl + C^P R )x c % + H.C., (19) 


*=ij=i 


such that 


C^J = gi-KrU&D^DZj - ^U* 2 D^ 3a D^ - n e U* 2 D t EAoi 
+ UiiD^DZj - K N U* 2 D% 2a D" 2j ) 


( 20 ) 


C«j =g{-K v V i 2 Dl\ ol b" j - K v V i 2 DZ a D^ - K„ e V i 2 DZ a D" 8j + V a Dl\ a D^ 
+ VnDlZD^ + V a DZ a D" 7j - K E V i 2 DZ 2 a D v 4j ), 


( 21 ) 


with 


(hn , H e ( 


(m.N, m T , m^, m e ) 
V2m w cos /3 


( 22 ) 


(he, 




.) = 


(m E ,m VT ,m I s fi ,m l , e ) 
V2mw sin f3 


(23) 
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and 


U*M C V = diag(m--,m--). (24) 

Next we discuss the interactions in the mass diagonal basis involving charged leptons, sleptons and 
neutralinos. Thus we have 
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X! H f ct{C^jP L + C' ( g j P R )Xif j + H.c., 
*=1 i =1 


such that 


where 


C'aij =V2{a Ti DZA - 5 Ei D™ a D T 2j - 7 riD^A + pEiD^A + 


- 7 nD^Dlj + aeiD^Dh - leiD^D^) 


C' a % APnDZAj - l Ei Dl* 2a D T 2j - 6 TZ Dl* la D T 3j + a Ei DZ 2 A + 


- ^niD E * 3a Dlj + 


PeiD E * 4a D^j — SeiDZA)’ 


grriEXlt 

®Ei = x-:—x , 

ztow sin p 

_ 5sin 2 6»rw* 

7-E« — e ^-li /a 5 

COS 


fe = e ^ + ^^ X2i (}" sin2(V 


$Ei = — 


gvi E X 4i 

2mw sin p ’ 


gm T X 3i 

°^ri — 0 o > 

2 mvF cos p 
§ = gmrX^ 

2 rriw cos P 


and where 


_ gm^X 3i 
a ^ 1 2m w cos (3 ’ 
c . = 9%^ 

2mu/ cos /3 ’ 


<5ei = - 


gm e X 3 i 
2mw cos P ’ 
gmeXZi 


2mw cos p ’ 


4ri — /7/ / ? Pei — 


COS #vy 


-^2* l^o+ sin ^ > 


Tti — 7m* — 7ei — 


, ssin 2 0 w , 

“CAi; i----AO? * 

* cos e w 


Here X' are defined by 


X' u = Xij cos Ow + X 2 i sin Q w 


XL = — Xu sin Ow + X 2 i cos O w ■ 



where X diagonalizes the neutralino mass matrix and is defined by the relation 


X T M^oX = diag(m x o,m x o,m x o,m x o) . (36) 

4 Interaction of leptons and mirrors with W and Z bosons 

In addition to the computation of the supersymmetric loop diagrams, we compute the contributions arising 
from the exchange of the W and Z bosons and the leptons and the mirror leptons in the loops. The relevant 
interactions needed are given below. For the W boson exchange the interactions that enter are given by 
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-C TW ^ = **££ ^ P \- C Z a p L + C% a Pn\r a + H.c., (37) 

i =1 a—1 

where 


r_ _9_\P)V* JJT | P)V* 7-)T , J-p* P)T 1 

° Lio ~~ \/2 + + ^Lli^LAoti > 

SIIV _ 9 \ T^iV* 7 ~)T 1 

° R ia — [ 1J R2i 1J R2a\ • 

For the Z boson exchange the interactions that enter are given by 


(38) 

(39) 


-£ 


ttZ 


= Z P £a=1 ELi T*Y[C?Pl + CiP R ]T P , 


where 




r>T 


dTLd 


La2 ±y L2/3 


DTLdT,* + D T 'D 


'La3^L30 


T t 


LotA 1 ^ LAfi / 


> 


,rf 

Lal^Ll/3 




DT La3 D L3P + 


1 T t 


and 



+ D Ra2 D R2P + D R a 3 D R3/3 + D 


yt 


4 r t 


os 9 


w 


l A«4 £) W) 


-r(d t 1.oD 


l Ra2 ly R2B) 


where x = sin 2 9w 


(40) 


(41) 


(42) 
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fi , 



Figure 1: The diagrams that allow the decay of p —> e + 7 via supersymmetric loops involving the chargino 
(top left) and the neutralino (top right) and via W loop (bottom left) and Z loop (bottom right) with 
emission of the photon from the charged particle inside the loop. 


5 The analysis of /i — >■ e + 7 Branching Ratio 


The decay p —> e + 7 is induced by one-loop electric and magnetic transition dipole moments, which arise 
from the diagrams of Fig(5] For an incoming muon of momentum p and a resulting electron of momentum 
p', we define the amplitude 

< e(p')\J a \p{p) >= u e (p')TaU^p ), (43) 


where 


r ^ = F^ e (q)ia afi q 0 + Fg e {q)(j a p 75 ^ + 


(44) 


m M +m e rrifj, + m e 

with q = p' — p and where m/ denotes the mass of the fermion /. The branching ratio of p —> e + 7 is given 

by 

24?r 2 


B{p -)■ e + 7 ) = 


G 2 F m 2 Jm M + TO e ) 


2{iir(o)r+iir(o)i 2 }, 


(45) 


where the form factors F% e and F£ e arise from the chargino, neutralino and vector bosons contributions as 


0 ) = ^; + + Fgo + F& + f% , 
*T ( 0) = F£ + + F£ 0 + F& + F£ . 


follows 


(46) 

(47) 














It is also useful to define B m and B e as follows 


Bra ( / ^ 

B e (n 


' + 7 ) = 


e + 7 ) = 


24t r 2 

GpTn^rrifi +m e ) 
24t r 2 

G 2 F m 2 (m F +m e ) 2 


r|^T e (0)| 2 , 


l*T(o)l 2 , 


(48) 

(49) 


where £> m is the branching ratio from the magnetic dipole operator and B e is the branching ratio from 
the electric dipole operator. We discuss now the individual contributions to F.F^ e supersymmetric and 
non-supersymmetric loops. 


The chargino contribution F^ + is given by 


2 8 r 

t= 1 3 =1 


m (.K + m e) f r.i /oL* , r-fl ^ M ’= 


' /✓''f.L /~iL* 1 rrti s~irC* \ tp 

‘ ^4i7 ^317 


647T 2 m 2 l v> 4ij'^3ij ~ ^4ij^Sij J 

Xi + V""; 

^e) r /^.R* \ 771 

64ir 2 m^.+ ^ ;i " + 1 m 2 . 


where 


and 


= 


1 


{—2a ; 3 — 3a ; 2 + 6 a; — 1 + 6 a : 2 In a;} 
F 2 (x) = -— {3a ; 2 — 4a; + 1 — 2a ; 2 lna;} . 


3(a: - l ) 4 

1 


(x-l)3 


The neutralino contribution F^ e 0 is given by 


2 .X 1 

4 8 


f 5=EE 

*=1 7=1 


-m F (m F + m e ) , L , Lt , R , R ( M l 


1927r 2 m 2 - 


f s~i!L /-yfL* 1 ri \ 771 

"t°4ij°3ij °4ij°3ij S r 3 


9 

ra ~ 0 
, Xi , 


( m A» + rn e) sn ,L n lR* , ,-V.R r" L *\p ( M- 
l°4ij L '3ii °4i i? °3ii 1-^4 


647T 2 TO^.o 


2 

ra ~ 0 

Xi° , 


where 


and 


F 3 (x) = 


(x~ 1 )‘ 


F 4 (x) = 


(*- l ) 3 

The contributions from the W exchange F.^y is given by 


-{—a : 3 + 6 a : 2 — 3a; — 2 — 6 a; In x} 


{—a ; 2 + 1 + 2 a; In a;} . 


y R e _ \ " )7 v( rr v + m e) r nW nW* , n W n 

r 2W - 327r 2 m^ ; ’ 1 mdL 


r/^vVy /tW* _i_ rrW s~iW*] tp I V’i 

[°Li4°Li3 °.Ri4 L '.Ri3J- r W r I o 

V , 


, 771^(771^ -h m e ) [ryW x 

+ 32TT*m* [G « 4 °« i3 + ^2 


(50) 

(51) 

(52) 


(53) 

(54) 

(55) 


(56) 
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where the form factors are given by 


and 


Fw(x ) = —-—r [4a ; 4 — 49a : 3 + 18a; 3 In a; + 78a; 2 — 43a; + 10] 

6 (x — l ) 4 


Gw(x) = -vr [4 — 15x + 12 x 2 — x 3 — 6 x 2 lnxl . 

(x — l ) 13 L J 


The contribution Fjfg from the Z exchange is given by 


where 


and 


we _ + m e) rrrZ n z* , n z nZ* 1 p ( 

*2Z - 647r 2 TO | [G Lf)4 C L03 + C R p 4 C Rfi3 \t z j 

, + m-e) lnZ r<z* , r^z n z* i/~< ( 

+ 647r 2 m| ^^ C R ^ + C R04 C L03 \G Z ^J , 

Fz(x) = ————vr [—5x 4 + 14x 3 — 39x 2 + 18x 2 lnx + 38x — 8] 

3{x — l) 4 


2 

Gz{x) = 7 - —77 lx 3 + 3x - 6 x In x — 4] 

(x — l ) 3 L J 

The chargino contribution F^ + is given by 


we 

3 X + 


where 


The neutralino contribution F^ 0 is given by 


V-V- ("V + m e )m^.+ (riL ^ „ R ^ ian ^ (m 2 - i+ ^ 

- 32 7 r 2 M? |C 4 ^ C 3 b - ( b, ( a-,!/ - ~^T » 

i=l 7 = 1 "1 \ v i / 

_ . . 1 r „ 2 lnx, 

fs(l) ' 2 (^TF { ' a: + 3 + rv7 1 - 


PF \ ' \ ' ( ?7 V + m e) m xi° tr"L r"R* r<’R c ( 

*3X° ~ Z^ Z^ 32tt 2 M 2 ' '-a 3,J ( W/ ( 3ij I *6 ^ 

i=l j=l T 1 \ 


M? 

Fj 


where 

r, / \ 1 r 2x In a; 

^e(a;) — XT-TZ^ l-} • 

2 (x — l) 2 1 — x 

The W boson contrition F 3 ^ is given by 


we _ _' 

1 QW - 


TO y>i ( m M + m e) r n w riW 


\/~lW /^ 1 V* S~iW S^W*1 T I ^ b j ] 
F'LiA'^Rtt L/ fli4°Li3 J 4 1 2“ I ’ 

\ m W / 


3 U Z-* 327 r 2 m.? v 

i—1 vv 


(57) 

(58) 


(59) 

(60) 

(61) 

(62) 

(63) 

(64) 

(65) 

( 66 ) 
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where the functions C} v and Cjj are given in section 4 and the form factor I\ is given by 


h{x) 


( I -*) 2 


( 11 1 2 

1- —X + -x 

4 4 


3a; 2 In x 
2(1 — x) 


And finally, the Z exchange diagram contribution F ^ is given by 


( 67 ) 


z?Me __ ( m M + me ) TOr ^ r/^-Z r<z* r>z nZ* 1 r I TOr a | 
3Z - Z. 327 T 2 rn | [ i4/3 ^3/5 C it4/3 C L3/3J7 2 ^ J ’ 


where the form factor I 2 is given by 

h(x) = 


(1 — a;) 2 

6 Estimate of size of B(fi —>■ e'j) 


1 1 o 3a; In x 

1 H —-x + —x + -j- -\ 

4 4 2(1 — x)_ 


( 68 ) 


(69) 


In this section we give a numerical analysis for the branching ratio £>(/i —► ey) . The analysis is done in 
an MSSM extension with soft breaking parameters taken at the electroweak scale. Thus no renormalization 
group running of GUT scale parameters is needed. The parameters entering the analysis are summarized in 
the appendix. The scalar mass and trilinear coupling parameters are mo and Aq in the slepton mass squared 
matrix. The corresponding ones in the sneutrino mass squared matrix are rrig and Aq where 

ml = M 2 tL = M 2 = Ml = Ml = M 2 L = Ml = M 2 eL = M e 2 
ml = M 2 n = M 2 Vt = Ml = Ml 
Aq = A t = A^ = A e = Ae 

Aq = A Vt = = A Ur = A n ( 70 ) 


The branching ratio B([i —> ey) arises as a consequence of mixing induced by the parameters fo, f” and 
fi, / 4 , fl where /’s are complex parameters and their arguments are the CP violating phases. The branching 
ratio —> ey) is a sensitive function of both the magnitudes as well as the phases of the mixing parameters 
/. We discuss the dependence of B(fi—> ey) on these below. 

Fig. ! exhibits the variation of £>(/r —► ey) as a function of the CP violating phases x '3 (left panel) and 
X 3 (right panel). As the two panels of Fig{2]show B(/i —> ey) is a sensitive function of these phases and 
can vary by an order of magnitude or more as the phases vary. The solid horizontal line gives the current 
experimental upper limit on B([i —> ey) from the MEG experiment [T. A very similar analysis holds when 
we vary the CP phases Xi an d y' 4 as exhibited in Fig{3j Figure [4] gives the relative strength of the magnetic 
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and the electric dipole transition operators to B(n —¥ ej) . Thus the left panel of Fig j4] exhibits the relative 
strength of the contributions from the magnetic dipole operator and the electric dipole moment operator 
as a function of the CP phase x!l- The right panel of Fig|4] exhibits the dependence of the electron EDM 
d e as a function of Xa where the solid horizontal line gives the current upper limit on d e from the ACME 
Collaboration (2HJ • Thus the right panel delineates the allowed regions of the parameter space, i.e., regions 
consistent with the experimental upper limit constraint on d e . The left panel of Fig[5] exhibits the variation 
of B(fi e'y) as a function of x!l f° r different values of the mixing parameter /" while the right panel of 
Fig[5] gives the electric dipole moment of the electron with the horizontal solid line giving the experimental 
upper limit on it. A comparison of the left and the right panels show the regions of x'i consistent with the 
current experimental upper limits on £>(/r —¥ e'y) and on d e and accessible with reasonable improvement in 
the sensitivity of experiment in the future. 

In tablet we illustrate numerically the relative strengths of the magnetic and the electric dipole transition 
operators to B(yi —> e'y) . Here we also show that the analysis is consistent with the current upper limits 
on the B(t —> ivy), on d e and the data on the neutrino masses. Thus the current experimental limit on 
B(t —> ny) is B{t —> n'y) < 4.4 x 10 -8 (BaBar) [2BJ and B(r —> ny) < 4.5 x 10 -8 (Belle) [IT]- Since the 
theoretical prediction of £>(r —> /jy) in this case is smaller by several orders of magnitude than the current 
experimental limit this decay mode is not of imminent interest in this case. In table 2 we give a numerical 
analysis of the form factors F% and F 3 and their sub pieces arising from the supersymmetric and the non- 
supersymmetric loops. Also listed are B m and B e as well as d e and the neutrino masses. One finds that 
typically the magnetic dipole contributions dominate the electric dipole contributions. The neutrino mass 
results of table 1 and table 2 are consistent with the constraint on the sum of the neutrino masses from 
cosmology, i.e., JAto^ < 0.44eV (95% CL) [25j and with the data on neutrino oscillations which give the 
neutrino mass squared differences so that [151 

A mh =ml-m\ = 2.4±°;^ x 1CT 3 eU 2 , (71) 

Anijj = ml - ml = 7.65^;^ x 1CT 5 eU 2 . (72) 

Fig. (|6| exhibits a variation of B(yi —» ey) as a function of the mirror masses m,E and m^. All points of 
the four curves of Fig. (6) are consistent with the constraints set on the neutrino masses by ?? and (??). 
We note that in the analysis of Fig. (2)- Fig.(4) the mass parameters are typically low. For instance in the 
analysis of Fig.4 we have used |^| = 310 GeV, Mi = 180 GeV and M 2 = 140 GeV and our vectorlike masses 
are chosen so that rag = tojv = 150,200,250,300. Such a choice may be close to the LHC exclusion plots 
based on LHC RUN I data and could be close to being probed with more data. We should note that the LHC 
particle searches are very model dependent as can be seen from the analyses of (31, 32 . For instance in the 
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ATLAS analysis of [32] the lightest chargino mass is excluded up to 700 GeV, 380 GeV, 345 GeV or 148 GeV 
for a massless neutralino depending on the allowed decay channels. These results would be even more model 
dependent if the neutralino is assumed massive with a varying mass. Thus while the current limits from 
LHC do not directly apply to our analysis, the choice of low mass parameters point to the possibility that 
they could be probed in RUN II of the LHC. It would thus be very interesting to carryout a signal analysis 
of this model specifically, for instance, for multilepton searches. Such an analysis, however, is beyond the 
scope of this paper. 

7 Conclusion 

In this work we have given an analysis of /i —> ey decay with inclusion of a vectorlike leptonic generation 
where mixings appear between leptons and mirror leptons as well as between sleptons and mirror sleptons. 
The decay /i —»• ey arises from diagrams with charginos and sneutrinos and mirror sneutrinos, and neu- 
tralinos, sleptons and mirror sleptons in the loops. Additionally electroweak loops are included where W, 
Z and leptons and mirror leptons and neutrinos are exchanged. An analytic analysis of these contributions 
is given in section 5 while a detailed numerical analysis is given in section 6 . Here it is shown that the 
current experimental limits from the MEG experiment put constraints on the parameter space of models. 
Further, the size of the new contributions are such that improvement in experiment will either reveal new 
physics or the improved experimental results will be able to probe large parts of the parameter space of the 
extended MSSM model. Thus the MEG experiment is continuing to collect data and is expected to explore 
the /x —> e + 7 decay down to a branching ratio sensitivity of a few times 10 “ 13 in the next few years. This 
will allow a further probe of this new class of MSSM extensions. 

Acknowledgments: PN’s research is supported in part by the NSF grant PHY-1314774. 


8 Appendix: The scalar mass squared matrices 

For convenience we collect here all the contributions to the scalar mass squared matrices arising from the F 
terms. They are given by 

^mass _ ^mass _|_ ^mass ( 73 ) 
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where £™ ass gives the mass terms for the charged sleptons while C^ ass gives the mass terms for the sneutrinos. 
For £™ ass we have 


/*mass 

~ L C 


r 1 + IM 2 + l/sl 2 + 1/3 I 2 ) E R E* r + ( + | /4 |2 + lyr/ja + v L El 


+ |/4|“ ) Tr? r + 


^1 l/l I 2 , I , | 2 A - , fvf\hl\ 2 , \ rf \2 


+ 1 / 3 1“ ft ft + 


1 / 4 1 2 f L R 9 *R 


vlM 2 


+ 1/3 1 AlAl + 


^ll^l 2 , , , //|2 \ vl\h 2 \ 2 , , ,„\ 2 


+ 1/4 1 2 Zr£r + 


+ 1/3? 


J /lM *«2 _ M*«2 _ f 2 H*Vi p p* / / 2 U 2 / 3 * fiVlf? \ p 

+ " -cr TtT * - - ^r £tE " + i^r + ^rj £tTi 


~7r + J £kt “ ++ El “ l+ 


V2 V2 j 1 
( fM* fkv iK 

V V 2 y /2 


ErFr 


fz*V2f' 2 , /M 


y/2 y/2 


-) ^ei + ( f -^- + ) E R e* R + f'f;a L r L + hf;^ 


+uft 


7 /* ~ ~ * , pH /»* ~ ~ * I pH p/* ~ . /»/ pH*~ ^2^2 ~ . 7-7 i 

4 e -R r i? + j3 /3 e £ r L + J3 J 3 e L^L + J4/4 e -RMi? — e -E/ e .R + H.C. > . 


We define the scalar mass squared matrix M~ in the basis (tl, El, tr, Er, JIl, ftR, e#)- We label the 
matrix elements of these as ( M~)ij = Mfj where the elements of the matrix are given by 

Ml 2 ! = M 2 L + + |/ 3 | 2 - m| cos 2/? Q - sin 2 <v) , 

V 22 = "+—yy—I" I/ 4 I 2 + I/ll 2 + I /4 | 2 + m l cos2/3sin“ 9 Wi 

M 33 = M 2 + + |/ 4 | 2 - m| cos2/3sin 2 0 W , 

M | 4 = M 2 + -yy— V I /3 1 2 + I /3 1 2 + I /3 |“ + m| cos 2/3 — sin 2 9 w ^j , 


M55 = M 2 l + —^ + I/ 3 1 2 - m| cos 2/3 Q - sin 2 9 w ^j , 
M | 6 = M 2 + ^-yt- + I/ 4 1 2 - m| cos 2/3 sin 2 0 W , 

M ^ 7 = M 2 l + + |/"| 2 - m| cos 2/3 Q - sin 2 9 w ^j , 

M 2 8 = M 2 + + |/ri 2 - m| cos2/3sin 2 0^ • 
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M | 2 = M|* = v2 ^ 3 


, ^13 = *4* = - »V 2 ), M 2 4 = M|* = 0, 


^lf 4 /1 ti^- 2 _ B/f2* _ fl 

VT~ T ~75T’ 13 “ 31 _ 72 

m 2 5 = M|* = /'/*, m 2 * = Mir = 0 , Mf; = m 2 ; = ttfiMl; = m|; = 0 , 


x/* 

J 2 


M | 3 = M 32 = 0, M | 4 = Mf = - m*i), M | 5 = m|; = 


V 2 / 3 / 2 * , *lfcl/4* 


\/2 


Mt = M? 9 * = 0, Mt = Mf = 


2 26 


*62 


27 


2 * 


* 72 


V2/3/2* 




, m 2 2 8 = m| 2 * = 0 , 


m 2 4 = m 4 2 3 * = 


U2/4/2 


v7 

m 3 2 7 = m 2 3 * = 0 , m 3 2 8 = m| 3 * = / 4 /r, 


\/2 

Vl/l/s *. m 3 2 5 = m 5 2 3 * = 0 , m | 6 = m 6 2 3 * = / 4 /r 


m | 5 = Mil = 0 , m | 6 = m 6 2 4 * = 


m 2 7 = Mt* = 0 , m/s = mH = 


*2/2/4* , *1/3^1 


*74 


48 


*84 


v7 

*2/2/4' 

\/2 


V 2 ’ 

* 1 / 3^2 




M 5 2 6 = M 6 t = ^=(* 1 ^ - M* 2 ), M 5 2 7 = Mt = t7r 


M 2 = Mi j? = 0, M 2 = Mil = 0 


J 58 

f 2 


*85 

f 2 * 


67 _ iV ^ 76 

2 _ !\/f 2 * 




m 6 2 8 = Mt = fjr,Mi s = Mi; = -^( Vl a; - »v 2 ). 


V 2 


We can diagonalize this hermitian mass squared matrix by the unitary transformation 


D t 'M?D t - *a 5 (M 2 M| 2 ,M| 3 ,M?,M| 5 ,M?,M| 7 ,M: 


/? ) 

T8' 


(75) 
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For C^ ass we have 


"TV 


Vl\f2? 


+ \fs\ 2 + \m 2 + \f^)N R N* R 

^r 1 + 1 / 5 1 2 + I/ 5 I 2 + l/s'l 2 ) N L N* L + ( V -M^ + |/ 5 | 2 ) v tR v* tR 


*>?l/(| 2 


V 2 o\h[\ 2 


2U1I 1 I r \2 \ ~ , I u 2 \' c l I , I - ~* , ( V 2.\h'l\~ 1 j,/, 2 \ ~ 

+ 1/3 1 ) VtLV tL + ( o 1 “ 1/3 1 ) + I - h I/5I J V^rV^r 


+ I + l/ 3 'l 2 ) wU + (77 + l/s I 2 1 KrKr 


f 2 V*V 2 £ T Ar * f[n*v 1_ _* A>*w 1. , (hvif'i /2U1/3* 


t ^2 * 
fhvift Kv 2 fS 


N l N* r - j -±^— VtLVtR 


72 " T "‘' TK 72 + V 72 72 

* , fW* &V l/n Ar (ftvrfZ f'i*V2h' 2 \ ~ 

NrUtR + - 7 TJ ^ + 1^71 - 7 l“j ' 


A7A 


V 72 


72 


, (h?V 2 n fs*V l/ 2 \ fr -* , / / 5 V 1/2 Ar 

+ 1.^1 - TO L eL + INI - J 


+ fsfP^rl + hfs*i>nRKR - 


h' 2 [JL*Vi 

-ppr' 1 


eL v e R 


+ fZf&eLKL + hfPVeRKR + fUs^L + fUP^R^R + #.C. j. 

Next we write the mass squared matrix in the sneutrino sector the basis (v t l,Nr, v t r, Nr, D^r, v^r, v e L , 7-r)- 
Thus here we denote the sneutrino mass squared matrix in the form ( M 2 )ij = m ?■ where 


"ill = 

II 

+ 

"1I7I 2 

2 

+ IM 2 

+ ^ m z cos 2/3, 

m 22 = 

- M 'n 

+ 

"II/ 2 I 2 

2 

+ l/ 5 | 2 - 
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+ I/ 5 F 

m 33 = 

= Ml 

+ 

vl\K | 2 
2 

+ IM 2 . 



^ to 

= m 2 x - 

,v 2 \f 2 \ 2 

2 

H- l/s| 2 -b l/sl 2 
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m 65 = 

= Kl 

+ 

"II7I 2 

2 

+ I7I 2 

1 

+ 2 "i 

■1 cos 2/3, 

m 66 = 

-K 

+ 

"II All 2 
2 

+ I7I 2 . 
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- MlL 

+ 

"IIAII 2 

2 

+ I/s I 2 

+ cos 2/3, 

00 to 

00 

= M 2 e 

+ 

"I a 2 | 2 
2 

+ 1/sT 

? 



16 



_2 


2* 


m 12 


m 2 1 


,2 


_2* 


m 14 


m 41 


2 


2* 


m 17 

— 

m 71 

— 

2 


_2* 


™23 


™32 


,2 


_2* 


™25 


™52 


2 


_2* 


m 26 


™62 


2 


,2* 


m 28 


™82 


2 


_2* 


m 35 


™53 


2 


2* 


m 37 


m 73 


m 45 

= 

2 * 

m 54 

= 

_2 


^,2* 


m 47 


m 74 


,2 


2 * 


™56 


™65 


,2 


^2* 


m 57 


m 75 


_2 


,2* 


m 67 


m 76 


,2 


,^2* 


m 78 

= 

m 87 

— 


V2 hf? 

V2 


^ 1 / 2/3 2 2 * fi ( a* \ 

^ > m i3 = m 3i = y|( V -* A -m), 


*■15 


2 * 

‘.51 


/343 ' 


‘16 


2 * 

'61 


2 2* /I 

. W 24 = ™42 = 

Vlftfs . KV2&* 

sft sft 

™2 _ „2* 


? 2 

'34 


? 2 

‘36 


*'43 


2* 

‘63 


2 2 * 

i 38 = m 83 

,2 _ 2 * 


,2 

M8 


7 2 * 

^84 


Vl f * 2/3 

, Kvzft 

72 

72 
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We can diagonalize the sneutrino mass square matrix by the unitary transformation 
D* MlD v = diag(M? , M? , M? , M? 4 , M? , M? , M? , M? ) . 
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Figure 2: Left Panel: An exhibition of the dependence of S(/x —> e'y) on x 3 where X 3 = X 3 = X 3 ■ The 
B(fi —> e'y) curves (bottom to top at X 3 = tt) are for the cases when |/ 3 | = |/ 3 | = \f^'\ = 5 x 10~ 5 , 7 x 10~ 5 , 
9 x 10 -5 , and 11 x 10~ 5 . Right Panel: An exhibition of B(yx —> ej) as a function of the X 3 where the 
curves (bottom to top at x 3 = tt) are for the cases |/ 3 | = 5 x 10~ 5 , 7 x 10 -5 , 9 x 10 —5 , 11 x 10 -5 where 
I/3I = I/3I = 5 x 10 -5 and %3 = X 3 = 0-3- The common parameters for both panels are tan/3 = 5, |/r| = 500, 
|Mi| = 130, |M 2 1 = 110, m N = 260, m E = 280, m 0 = 4x 10 5 , mg = 5x 10 5 , |A 0 | = |Ag| = 6x 10 5 , = 0.4, 
a 2 = a^ = 0.5, cxa 0 = a A > = 1, |/ 4 | = \fi\ = 1, \f(\ = 0.1, |/ 5 | = 3 x 10" 6 , |/^| = 8 x 10" 7 , |/^| = 5 x 10~ 6 , 
X 4 = 1, X '4 = X 4 = 0-5, X 5 = X 5 = X 5 = 1- The solid horizontal line is the upper limit from the MEG 
experiment [T. Here and in the rest of the figures and in the tables all masses are in GeV and phase angles 
in radian. 
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Figure 3: Left Panel: Plot of B(n e'y) as a function of X 4 when x 4 = X 4 = X 4 - The curves from bottom 
to top at X 4 = 7T are for the cases |/ 4 | = I/ 4 I = I /4 | = 0.1 , 0.5, 1, and 1.5 and |/ 3 | = |/g| = |/g | = 5 x 10 -5 
and X 3 = X 3 = X 3 = 0-3- The solid horizontal line is the upper limit from the MEG experiment JT]. Right 
Panel: An exhibition of the dependence of B(n —» e'y) on x!a- The curves from bottom to top at X4 = 7r are 
for the cases when |I =0.3, 0.4, 0.5, and 0.6 and |/ 4 | = |/"| = 0.4,x 4 = X 4 = 1- All other parameters in 
both panels are the same as in Fig(2] 
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Figure 4: Left Panel: An exhibition of the magnetic contribution B m (dotted) as given by Eq. (??), the 
electric contribution B e (dashed) as given by Eq. (??) to B(^l —► ey) where the solid curve stands for their 
sum as a function of x 4 when title = m/v = 300. Right Panel: An exhibition of d e as a function of x 4 • 
The curves from top to bottom at y 4 = tt are for values of to# = tun = 150, 200, 250, and 300. The 
solid horizontal line is the upper limit on d e from the ACME Collaboration p[5j. The common parameters 
for the two panels are |/t| = 310, \M\\ = 180, |M 2 | = 140, tan/3 = 20, m 0 = 4 x 10 5 , |Ao| = 1.5 x 10 6 , 
ttiq = 4x 10 5 , |Aq| = 5.1 x 10 6 , ay = 0.4, a 2 = 0.2, = 0.7, aA 0 = Ua* = 1. The mixings are |/ 3 | = 7x 10” 4 , 

|/'| = lxl0- 4 , |/"| =2xl0- 4 , |/ 4 | =3xl0- 2 , |/1| = 0.4, |/"| = 5xl0- 2 , |/ 5 | = 3.8x HT 6 , |/'| = 2.2xl0” 6 , 
i/s I = 3 x 10“ 6 , X 3 = X 3 = X3 = 1, X4 = 0.3, x'4 = 0.2, xs = X5 = X5 = 0-5. 



(i) rriE = tun = 300 

(ii) rriE = rriN = 150 

Bm 

1.1 x 10~ 12 

8.3 x 10~ 12 

B e 

1.6 x 10~ 13 

1.2 x 10 -12 

B{/J- ey) 

1.2 x 10~ 12 

9.4 x 10” 12 

B(t /xy) 

4.8 x 10 -25 

8.0 x 10- 25 

d e (ecm) 

6.3 x 10" 30 

1.3 x 10 -29 

m v 3 

5.0 x 10” 11 

5.5 x KT 11 

m u2 

8.9 x 10” 12 

9.6 x KT 12 

m v 1 

1.1 x 10~ 12 

2.5 x 10~ 12 


Table 1: An exhibition of the numerical values of B rn , B e , B(n ey) and B{t — > /z + 7 ) when y 4 = 2 for 
two values of = tun while other parameters are the same as in Fig[4] The values of the electron EDM 
d e and the neutrino masses m Ul , m„ 2 , m„ 3 are also exhibited. 
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Figure 5: Left Panel: A plot of B(/j, —l ex) as a function of X3 where the curves from bottom to top at x" = 2 
are for values of 1/3'| = 7 x 10” 5 , 8 x 10~ 5 , 9 x 10 -5 , 1 x 10 -4 . The solid horizontal line is the upper limit from 
the MEG experiment pQ. Right Panel: Plot of d e as a function of x '3 for the same values of /" as in the left 
panel. The solid horizontal line is the upper limit on d e from the ACME Collaboration [23]. The parameters 
used are tan/? = 5, \fi\ = 500, || = 130 , |M 2 | = 110 , TOjv = 260 , rriE = 240 , m 0 = 7 x 10 4 , ijIq = 5 x 
10 4 , |Ao| = \Aq\ = 6 x 10 5 , op = 0.4 , a 2 = = 0.5, o.a q = ol a » = 1. The mixings are I/3I = 3 x 

10- 5 , \fi\ = 4 x 10- 6 , |/ 4 | = I/ll = 0.8 , |/i'| = 0.1 , |/ 5 | = 3 x° 10- 6 , I/ll = 7 x 10~ 6 , |/ 5 "| = 5 x 
10 -6 . Their CP phases are %3 = 0.3 , X 3 ' = 0.4 , x ' 4 = 1 > X 4 = Xa = 0.5 - X5 = X5 = Xh" = 1- 


(i) X '3 = 0-39 (>i) \'(-2 1 

b 7 V 1 9 si V 1 fpTBT=XTs 


pfie 

4 2 X + 

F fie 

r 2x a 

t 2W 
^2 Z 

2.7 x 10~ i8 e _u - Y4 * 

3.4 x 10 _21 e°- 56i 

5.8 x 10 _15 e _1 ' 5i 

4.0 x 10 15 e 1 - 4i 

2.8 x 10 _iS e _2 i * 

3.5 x 10- 21 e-°- 391 

6.1 x 10 _15 e ol2i 

4.7 x 10 15 e o l5i 

FH.0) 

1.9 x 10 -15 e -1 - 11 

1.1 x 10 14 e UUU24i 

Bm 

5.3 x 10~ 14 

5.6 x 10 -13 


2.4 x 10~ 18 e 3U4 “ 

9.5 x 10 19 e U 45 * 

F n e 

r 3x° 

1.3 x 10 _2O e°- 53i 

1.6 x 10 _21 e _l li 

1 3W 

3.3 x lO-^e- 1 - 71 

2.3 x io- 15 e -° 0191 

wxe 

^3Z 

1.5 x 10 15 e 1 - 9i 

5.8 x 10 16 e O 34i 

*T( 0) 

2.2 x KT 15 e~ 2 11 

2.9 x 10 _15 e O O51i 

B e 

6.6 x KT 14 

6.7 x 10“ 14 

->■ e l) 

1.2 x KT 13 

1.7 x 10 -12 

B(t iij) 

2.5 x KT 19 

2.5 x 10" 19 

d e (ecm) 

4.7 x 10~ 29 

6.1 x 10 3U 

m v 3 

4.7 x KT U 

4.7 x 10" 11 

m v 2 

9.1 x KT 12 

8.8 x 10~ 12 

m v 1 

2.1 x KT 12 

4.9 x 10 -13 


Table 2: An exhibition of the numerical values of the form factors F£ e and F^ e and their sub pieces for two 
cases: (i) and (ii). For case (i) |/g| = 10 -4 and x'3 = 0.39 while for case (ii) |/g | = 0.7 x 10 -4 and x'3 = 2.1. 
All other parameters used in this table are the same as Fig[5j The magnetic and electric transition operators 
B rn and B e are also listed as are d e and and the neutrino mass eigenstates for each case listed above. 
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Figure 6 : A plot of Z?(/i —► e'y) as a function of the = tun where the curves from top to bottom at 
rriE = niN = 150 are for values of I/ 4 1 = 0.1, 0.5, 1, 1.5. The parameters used are |/z| = 500, \M\\ = 130, 
\M 2 \ = 110, tan/3 = 10, m 0 = 4 x 10 4 , |A 0 | = 6 x 10 5 , mg = 5 x 10 4 , |Ag| = 6 x 10 5 , a x = 0.4, a 2 = 0.5, 
= 0.5, a Ao = a A p = 0.6. The mixings are |/ 3 | = 5 x 10" 5 , |/ 3 | = 5 x 10" 5 , |/ 3 | = 5 x 10" 5 , I/ 4 I = 0.5, 
\f’l\ = 5 x 10- 1 , \M = 3 x 10- 6 , \f 5 \ = 8 x 10- 7 , \fH\ = 5 x 10" 6 , X 3 = xg = X 3 = 0.3, X 4 = X 4 = X? = 1, 
X5 = X 5 = X5 = 1- 
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